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We  operated  our  coherent  microwave  radars  at  sea  a  total  of  four  times  in  this  project.  We 
operated  a  radar  we  call  RiverRad  on  three  ships,  the  RW  Revelle  and  RW  Ocean  Researcher  1 
in  the  South  China  Sea,  and  the  RW  Endeavor  off  the  coast  of  New  Jersey.  In  addition,  we 
operated  a  similar  airborne  radar,  CORAR,  on  a  Cessna  Skyhawk  off  the  New  Jersey  coast  at  the 
same  time  as  the  Endeavor  cruise.  A  subset  of  the  results  from  these  cruises  and  flights  has  been 
submitted  to  the  Journal  of  Geophysical  Research-Oceans  for  publication.  They  have  the 
following  titles  and  are  attached  to  this  report: 

“The  Normalized  Radar  Cross  Section  of  the  Sea  for  Backscatter:  Part  1 .  Mean  Eevels” 

“The  Normalized  Radar  Cross  Section  of  the  Sea  for  Backscatter:  Part  2.  Modulation  by  Internal 
Waves” 

In  addition  to  these  two  papers,  a  third  paper  dealing  with  our  radar  measurements  of  internal 
wave  characteristics  in  the  South  China  Sea  is  in  preparation  and  will  be  submitted  shortly.  Its 
title  is  the  following: 

“Characteristics  of  Internal  Waves  in  the  South  China  Sea  Observed  by  Shipboard  Radar” 

Some  of  the  results  presented  in  these  papers  can  be  briefly  summarized  as  follows: 

•  The  multiscale  model  we  developed  previously  to  explain  the  dependence  of  the  normalized 
radar  cross  section  (NRCS)  of  the  sea  at  microwave  frequencies  on  wind  speed  and  direction  is 
viable  at  VV  polarization  at  much  larger  incidence  angles  than  we  had  thought.  At  this 
polarization  it  works  out  to  89  degree  incidence  angles.  By  contrast  at  HH  polarization  the 
model  underpredicts  the  NRCS  of  the  sea  for  incidence  angles  above  about  45  degrees. 

•  At  high  grazing  angles,  HH  polarized  cross  sections  maximize  upwind  and  minimize 
downwind;  upwind  they  are  slightly  smaller  than  VV  cross  sections  .  By  contrast,  at  VV 
polarization  cross  sections  maximize  upwind,  minimize  cross  wind,  and  show  a  second  local 
maximum  downwind,  a  more  standard  behavior. 


•  Generally  the  NRCS  for  W  exceeds  HH  as  predicted  by  theory.  However,  situations  exist  on 
the  ocean  where  HH  exceeds  VV  by  as  much  at  15  dB.  These  remain  unexplained. 
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•  Internal  waves  (IW)  modulate  the  mean  NRCS  to  produce  maximum  cross  sections  near  the 
IW  peak.  However  the  exact  location  of  the  maximum  varies  depending  on  the  strain  rate  of  the 
IW.  For  the  low  strain  rates  found  in  deep  water,  the  maximum  occurs  at  the  crest  while  for  the 
larger  strain  rates  often  found  in  shallow  water,  the  maximum  occurs  well  ahead  of  the  crest, 
sometimes  as  much  as  a  kilometer. 

•  Internal  waves  often  cause  the  HH  NRCS  to  maximize  at  larger  values  than  the  VV  NRCS, 
probably  due  to  breaking  wave  effects.  This  calls  into  question  the  historical  explanation  of  the 
effect  of  IWs  on  the  NRCS,  which  was  based  on  the  modulation  of  wind  waves.  Breaking  of 
short  gravity  waves  also  appears  to  be  involved. 

•  The  visibility  of  internal  waves  in  radar  imagery  appears  to  vary  little  with  azimuth  angle. 

•  In  addition  to  the  large,  non-linear  solitons  that  exist  in  the  South  China  Sea,  smaller  nearly 
sinusoidal  wave  trains  also  exist  that  propagate  at  speeds  different  than  those  of  the  solitons, 
interact  with  the  solitons,  and  are  modulated  by  the  solitons. 

•  The  solitons  travel  at  speeds  that  depend  on  their  amplitudes  in  a  manner  explained  very  well 
by  the  standard  two  layer  model.  These  speeds  reach  3  m/s  in  deep  water  and  just  over  1  m/s  on 
the  shelf 
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The  normalized  radar  cross  section  of  the  sea  for  backscatter,  NRCS  or  Oq,  is 
investigated  for  incidence  angles  between  0°  and  89°  using  a  variety  of  data  collected 
over  more  than  two  decades.  For  all  of  these  data,  the  radar  return  has  been  compared  to 
that  from  comer  reflectors  to  yield  calibrated  normalized  radar  cross  sections.  The  most 
recent  measurements  were  made  from  several  ships  using  a  coherent,  dual-polarized,  X- 
band  radar.  These  measurements  show  that  Oo  (VV)  at  high  incidence  angles  exhibits 
wind  speed  and  azimuth  angle  dependence  similar  to  those  at  the  lower  incidence  angles 
used  in  scatterometry.  It  is  nearly  as  large  looking  downwind  as  looking  upwind  and 
minimizes  near  the  crosswind  direction.  ao(HH)  behaves  very  differently.  At  incidence 
angles  of  88°  to  89°,  it  is  largest  looking  upwind  and  smallest  looking  downwind. 

Fits  of  the  multiscale  model  of  microwave  backscatter  from  the  ocean  (Plant,  2002)  to 
these  data  along  with  a  variety  of  data  collected  previously  at  lower  incidence  angles 
shows  that  over  the  whole  range  of  incidence  angles  from  0°  to  89°,  ao(W)  is  rather  well 
explained  by  the  model  while  measured  ao(HH)  values  are  generally  higher  than  the 
model  predicts  at  incidence  angles  above  about  45°.  This  suggests  that  scattering 
phenomena  exist  on  the  ocean  surface  that  affect  HH  backscatter  very  strongly  at  the 
higher  incidence  angles  while  impacting  VV  polarized  backscatter  only  slightly. 

This  suggestion  is  strengthened  by  our  observation  of  high- incidence-angle 
backscatter  from  the  ocean  where  ao(HH)  exceeds  ao(VV)  by  as  much  as  15  dB  for  long 
periods  of  time  and  over  large  spatial  areas.  We  examine  phenomena  that  might  account 
for  this  behavior,  including  bound,  tilted  Bragg  scatterers,  specular  points,  wedges,  spray, 
and  multipath  dihedral  and  cylinderical  features.  We  suggest  that  multipath  dihedral-type 
features  are  likely  to  be  relevant  scatterers  since  they  produce  large  ao(HH)/ao(VV)  due 
to  Brewster  damping  of  the  first  W  bounce. 

Introduction 

While  the  normalized  radar  cross  section  of  the  sea,  NRCS  or  Oo,  for  cross  polarization 
can  be  illuminating,  for  instance  in  detecting  multiple  scattering  events,  it  is  generally 
much  smaller  than  copolarized  cross  sections  (Wiltse  et  ah,  1957;  Kalmykow  and 
Pustovoytenko,  1976;  Lee  et  al,  1999).  For  this  reason,  Oo  values  obtained  with  the 
electric  field  of  the  radiation  either  vertical  ao(VV)  or  horizontal  ao(HH)  on  both 
transmission  and  reception  are  more  useful  modes  and  ones  that  have  been  studied  much 
more  than  cross  polarized  cross  sections.  In  this  series  of  two  papers  we  report  our  work 
on  the  copolarized  normalized  radar  cross  section  of  the  sea  for  backscatter,  Oq.  The 
present  paper.  Part  1 ,  concentrates  on  the  behavior  of  the  mean  Oo  as  a  function  of  wind 
speed,  incidence  angle,  and  polarization.  Part  2  will  document  our  recent  measurements 
on  Oo  modulated  by  surface  currents  set  up  by  internal  waves  in  the  ocean. 


Historically  ao(VV)  has  been  known  to  be  rather  well  explained  by  Bragg  seattering 
augmented  by  a  eomposite,  or  two-seale,  sea  surfaee  over  the  range  of  ineidenee  angles 
from  approximately  20°  to  60°  and  probably  at  even  higher  ineidenee  angles,  into  the  so- 
called  low-grazing  angle  regime  (Wright,  1968;  Bass  et  ah,  1968).  ao(HH),  on  the  other 
hand,  has  been,  and  eontinues  to  be,  more  mysterious.  It  has  appeared  to  be  fairly  well 
predieted  by  Bragg/eomposite-surfaee  seattering  from  20°  to  45°  ineidenee  but  exeeeds 
the  predietions  at  larger  ineidenee  angles.  At  lower  ineidenee  angles,  speeular  seattering 
has  been  used  to  aeeount  for  both  ao(VV)  and  ao(HH),  but  requires  an  “effective” 
reflection  coefficient.  Matching  specular  and  Bragg  scattering  near  20°  ineidenee  is  a 
non-trivial  exereise. 

In  2002,  Plant  developed  a  multiscale  model  and  showed  that  by  adding  a  third  seale 
of  the  sea  surfaee  and  using  a  Kirehoff  integral  approaeh,  whieh  is  well  approximated  by 
Bragg  seattering  for  low  wind  speeds  or  moderate  to  large  ineidenee  angles,  predietions 
in  agreement  with  ao(VV)  and  ao(HH)  measurements  eould  be  produeed  over  the 
ineidenee  angle  range  0°  to  about  45°  without  invoking  an  arbitrary  effeetive  reflection 
coefficient  (Plant,  2002).  Applying  the  theory  beyond  this  incidence  angle  range  was 
diseouraged  in  that  paper,  however,  beeause  multiple  scattering  and  breaking  waves  are 
not  ineluded  in  the  model. 

In  this  paper  we  show  that  this  multiseale  model  in  faet  explains  ao(VV)  over  the 
entire  range  of  ineidenee  angles  from  0°  to  89°  if  shadowing  is  ineluded.  ao(HH)  remains 
unexplained  by  the  model  beyond  an  ineidenee  angle  of  about  45°,  however.  We  will 
investigate  possible  reasons  for  this  behavior  of  ao(HH),  ineluding  baekseattering  from 
spray  and  from  surfaee  features  eaused  by  breaking  waves. 

Low-Grazing  Angle  Backscatter  from  the  Ocean 

An  X-band  Doppler  radar  was  operated  on  the  R/V  Revelle  in  the  South  China  Sea  in 
2005,  on  the  RW  Endeavor  off  the  New  Jersey  eoast  in  2006,  on  the  RW  Melville  in  the 
Philippine  Sea  in  2007.  The  radar  transmitted  a  peak  power  of  10  watts  and  used  a  pulse 
width  of  15  m  in  2005,  3.75  m  in  2006,  and  30  m  in  2007.  The  radar  is  the  RiverRad 
system  that  was  developed  for  river  surfaee  veloeity  measurements  and  is  deseribed  more 
fully  in  Plant  et  ah,  2005.  It  is  shown  mounted  on  the  three  ships  in  Figure  1 . 

In  all  three  years,  we  ealibrated  the  sea  return  against  that  from  a  eorner  refleetor  so 
that  fully  calibrated  normalized  radar  eross  seetions,  Oo,  eould  be  obtained.  In  2005,  only 
VV  polarized  baekseatter  was  obtained  due  to  a  switeh  failure  and  the  antennas  were 
fixed  looking  nearly  toward  the  bow  of  the  ship.  In  2006,  both  HH  and  VV  baekseatter 
were  eolleeted.  The  antennas  seanned  through  approximately  80°  and  were  directed  35° 
apart  in  azimuth.  In  2007,  the  antennas  were  again  fixed  but  this  time  looking 
perpendieular  to  the  heading  of  the  ship;  both  HH  and  VV  baekseatter  were  eolleeted. 

On  all  cruises,  baekseatter  was  eolleeted  both  from  seas  disturbed  only  by  the  wind  and 
from  seas  where  internal  waves  were  present. 


Figure  1.  RiverRad  mounted  on  the  ships,  a)  RA^  Revelle  in  2005,  b)  RA^  Endeavor  in 
2006,  c)  R/V  Melville  in  2007. 

Figure  2  shows  the  angular  dependence  of  ao(VV)  observed  in  2005  at  an  89° 
incidence  angle  and  at  various  wind  speeds  in  the  absence  of  internal  waves.  Note  that 
ao(VV)  exhibits  the  same  second-harmonic  dependence  on  azimuth  angle  that  it  does  at 
lower  incidence  angles.  This  is  not  the  case  for  the  horizontally  polarized  cross  section, 
ao(HH)  shown  in  Figure  3  at  nearly  the  same  incidence  angle.  Here  both  ao(VV)  and 
ao(HH)  are  shown  as  determined  from  the  measurements  of  2006.  ao(HH)  maximizes 
when  the  antenna  looks  into  the  wind  and  minimizes  in  the  opposite  direction.  These 
results  agree  with  the  uncalibrated,  nearshore  measurements  of  Trizna  and  Carlson,  1996. 
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Figure  2.  Angular  dependence  of  ao(VV)  at  X-band  and  an  incidence  angle  of  89° 

0.25°.  These  2005  results  are  for  the  ocean  surface  disturbed  only  by  wind.  Symbols  are 
for  various  wind  speeds;  +  =  4  m/s,  o  =  6  m/s,  +  =  8  m/s. 
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Figure  3.  a.  Angular  dependence  of  ao(W).  b.  Angular  dependence  of  ao(HH).  These 
2006  results  are  for  a  sea  surface  disturbed  only  by  wind;  the  incidence  angle  is  88.3° 
0.1°.  Symbols  as  in  Figure  2. 

Figures  4  and  5  show  the  wind  speed  dependences  of  these  data  and  compare  them  to 
simple  Bragg  scattering  predictions  using  the  spectrum  of  Elfouhaily  et  ah,  1997.  Clearly 
the  Bragg  predictions  are  much  too  low  and  the  increase  with  wind  speed  is  too  small.  In 
fact,  Bragg  predictions  for  ao(HH)  are  below  the  bottom  of  the  plot  in  Figure  5b. 


Wind  Speed,  m/s 

Figure  4.  Wind  speed  dependence  of  ao(VV)  at  various  azimuth  angles  from  the  wind  for 
the  data  of  2005.  y  is  the  angle  between  the  antenna  look  direction  and  the  direction  from 
which  the  wind  comes.  Lines  show  simple  Bragg  scattering  predictions  using  the  wind 
wave  spectrum  of  Elfouhaily  et  ah,  1997. 


Figure  5.  Wind  speed  dependenee  of  Oo  at  various  azimuth  angles  from  the  wind  for  the 
data  of  2006.  a.  ao(W),  b.  ao(HH).  Lines  show  simple  Bragg  seattering  predietions  for 
Oo(VV)  using  the  wind  wave  speetrum  of  Elfouhaily  et  ah,  1997;  ao(HH)  predieted  by 
Bragg  seattering  is  below  the  bottom  of  the  lower  graph. 

This  under  predietion  ean  be  overeome  for  ao(VY)  by  resorting  to  the  multiseale 
model  (Plant  2002).  This  model  extended  the  standard  eomposite  surfaee  model  to 
inelude  effeets  of  waves  that  were  intermediate  between  Bragg  seattering  waves  and 
long,  modulating  waves.  This  model  used  results  of  the  small  slope  approximation 
(Voronovieh,  1985)  and  the  integral  expansion  method  (Fung  et  al,  1992)  for  baekseatter 
to  eombine  speeular  and  Bragg  seattering  into  a  single  Kirehhoff  integral  that  eovered 
ineidenee  angles  from  nadir  to  mid-range.  The  predietions  of  this  model  were  shown  to 
agree  well  with  the  Ku-band  airship  data  mentioned  above  and  with  other  data  sets  at  C 
and  Ka  bands  as  well  as  Ku-band  at  nadir  ineidenee.  These  eomparisons,  however,  were 
only  done  for  ineidenee  angles  up  to  50°  beeause  breaking  waves  and  spray  were  not 
ineluded  in  the  model,  making  it  questionable  above  this  ineidenee  angle. 

The  improved  fit  using  this  model  is  demonstrated  in  Figure  6  where  the  left  panels 
show  ao(W)  versus  wind  speed  for  upwind,  downwind,  and  erosswind  looks.  The 
squares  in  this  figure  are  results  of  the  multiseale  model  modified  to  inelude  simple 
geometrie  shadowing.  Intermediate-seale  faeets  were  not  ineluded  in  the  ealeulations  if 
they  eould  not  be  seen  by  the  antenna.  Bragg  seattering  predietions  ean  also  fit  the  data  if 
a  mean  tilt  of  the  illuminated  parts  of  the  surfaee  given  by  U/4  is  assumed,  where  tilt  is  in 
degrees  and  U  is  in  m/s.  This  yields  the  red  lines  in  Figure  6 

Figure  6  also  shows  similar  data  and  models  for  ao(HH)  in  the  right  panels.  Here  the 
measured  values  are  mueh  larger  than  either  those  predieted  by  the  multiseale  model  or 
by  tilted  Bragg.  This  is  true  even  when  the  antenna  is  direeted  down  wind  and  wave. 
Sinee  breaking  wave  effeets  are  generally  assumed  to  have  little  effeet  on  baekseatter 


with  the  antenna  direeted  nearly  down  wave,  we  interpret  these  results  to  mean  that  a 
seattering  mechanism  exists  at  these  grazing  angles  other  than  multiple  scattering  from 
breaking  waves  that  affects  ao(HH)  much  more  than  ao(W).  This  is  reasonable 
considering  the  much  higher  level  of  ao(VV)  shown  in  the  downwind  part  of  Figures  3 
and  5.  Thus  we  need  an  increase  in  the  ao(HH)  values  predicted  by  the  multiscale  model 
of  about  10  dB  in  the  downwind  direction  and  about  20  dB  in  the  upwind  direction  in 
order  to  explain  the  data. 
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Figure  6.  Fit  of  the  multiscale  model  (squares)  to  the  data  of  2006  for  upwind  looks 
(circles).  The  red  lines  are  Bragg  scattering  predictions  if  the  local  grazing  angle  in 
degrees  is  larger  than  the  nominal  one  by  U/4,  where  U  is  wind  speed  in  m/s. 

Fits  to  the  Multiscale  Model  over  all  Incidence  Angles 

To  check  the  predictions  of  the  multiscale  model  further,  we  compared  it  with  lower 
incidence  data  in  addition  to  these  high-incidence-angle  measurements.  Over  the  last  two 
decades,  we  have  collected  a  variety  of  ao(VV)  and  ao(FlH)  data  at  14.0  GFlz  (Ku-band) 
and  9.36  GFlz  (X-band).  These  data  were  collected  at  many  different  incidence  angles, 
wind  speeds,  and  azimuth  angles  from  platforms  as  varied  as  research  platforms,  airships, 
airplanes,  and  ships  (Plant,  1997;  Plant  et  ah,  1998,2005). 

Figures  7-10  show  the  comparison  out  to  higher  incidence  angles  that  shown  in  Plant 
(2002).  The  predictions  of  the  multiscale  model  for  ao(VV)  and  ao(FlFl)  are  compared 
with  the  following  data  sets:  Ku-band  collected  on  the  German  Research  Platform 
Nordsee  (Plant  1997,  2003b),  on  the  airship  off  the  coast  of  Oregon  (Plant  et  ah,  1998), 
and  the  TOPEX/Poseidon  satellite  altimeter  (Flwang,  1998);  X-band  data  from  a  Twin 
Otter  aircraft  flown  off  the  coast  of  North  Carolina  (Plant  et  al.,2005)  and  the  ship  data 
discussed  above.  We  ran  the  model  with  various  ocean  wave  spectra  and  for  both  X  and 
Ku  band.  Figures  7-10  show  the  results  for  the  spectrum  that  Plant  (2002)  calls  the  D- 
spectrum,  a  combination  of  spectra  suggested  by  Donelan,  Banner,  and  Plant. 
Comparisons  of  predictions  with  other  spectral  models  can  change  the  model  predictions 
somewhat. 


Figures  7-10  show  predictions  of  the  multiscale  model  at  both  Ku  and  X-bands.  They 
are  very  similar  above  the  lowest  wind  speed,  differing  by  3  dB  at  most  for  a  4  m/s  wind 
speed.  The  data  at  wind  speeds  of  8  (Figure  8),  12  (Figure  9)  and  16  m/s  (Figure  10)  all 
indicate  that  the  model  predicts  ao(W)  rather  well  over  the  whole  range  of  incidence 
angles  from  0°  to  89“  for  upwind,  crosswind,  and  downwind  antenna  look  directions.  At 
4  m/s  (Figure  7),  the  fit  is  less  good  in  places  but  wave  spectra  at  this  wind  speed  can 
vary  widely  due  to  wind  variability  (Plant,  2000).  The  fit  is  not  as  good  for  ao(FIFI).  For 
this  polarization,  the  predictions  appear  to  be  accurate  up  to  an  incidence  angle  of  about 
45°  but  fall  significantly  below  the  data  for  higher  incidence  angles.  This  is  true  for  all 
three  antenna  look  directions.  In  a  following  section,  we  will  discuss  possible  reasons 
that  the  multiscale  model  might  under  predict  ao(HH)  at  high  incidence  angles. 


Figure  7.  Comparison  of  predictions  of  the  multiscale  model  with  data  for  a  wind  speed 
of  4  m/s.  The  D-Spectrum  of  Plant  (2002)  was  used  in  the  calculation.  Black  lines  are 
predictions  at  Ku-band  while  dashed  lines  are  predictions  at  X-band.  Red  always 
indicates  ao(HH)  while  black  is  ao(VY).  Symbols  for  data  are  as  follows:  circles  -  Ku- 
band  data  from  the  TOPEX/Poseiden  altimeter;  diamonds  -  Ku-band  from  the  airship 
(Plant  et  ah,  1998);  down  triangles  -  X-band  from  the  airplane  (Plant  et  ah,  2005);  up 
triangles  -  X-band  data  from  the  ships  discussed  above. 


Figure  8.  Same  as  Figure  7  but  for  a  wind  speed  of  8  m/s.  The  squares  are  Ku-band  data 
from  the  Nordsee  platform  (Plant,  1997,  2003b) 


Figure  10.  Same  as  Figures  7  and  8  but  for  a  wind  speed  of  16  m/s. 

Anomalous  Low-Grazing  Angle  Backscatter  from  the  Ocean 

In  Part  2  of  this  paper,  we  will  show  that  internal  waves  can  sometimes  generate 
surface  conditions  for  which  ao(HH)>  ao(W),  which  as  shown  above  does  not  normally 
occur  for  ocean  backscatter  on  large  scales.  However,  we  have  observed  an  occurrence 
of  mean  ao(HH)  larger  than  mean  ao(VV)  in  2007  when  the  RiverRad  operated  from  the 
R/V  Melville  around  the  Philippine  Islands.  In  the  images  collected,  we  observed  many 
regions  of  high  oceanic  backscatter  of  unknown  cause.  An  example  is  shown  in  Figure 
1 1 .  A  total  of  1 5  images  of  similar  quality  were  collected.  The  figure  shows  ao(HH)  and 
ao(VV)  in  dB  with  associated  scatterer  velocities  in  m/s.  Also  shown  in  the  bottom  panel 
of  the  figure  are  the  wind  speed,  wind  direction,  ship  speed,  and  ship  heading.  The 
antennas  looked  broadside  to  the  ship  on  the  starboard  side  and  were  16.5  m  above  the 
water.  The  high  cross  section  and  velocity  features  cover  more  than  2  km  in  range  and 
lasts  for  nearly  5  minutes  as  the  ship  goes  by  at  6  m/s.  Thus  the  features  are  also  about  2 
km  in  the  azimuth  direction.  ao(HH)  is  5  to  15  dB  larger  than  ao(VV)  over  most  of  this 
region. 

The  large  regions  covered  by  this  anomalous  backscatter  lead  us  to  believe  that  it  is 
not  of  man-made  origin,  for  instance  interference  from  other  marine  radars.  This  belief  is 
strengthened  by  the  fact  that  most  of  the  high-backscatter  regions  occur  near  the  ship  at 
the  same  time  that  shifts  in  wind  speed  or  direction  were  measured  by  the  ship’s 
anemometer.  We  cannot  rule  out  the  involvement  of  rainfall  in  the  occurrence  of  these 
features  since  precipitation  measurements  on  the  ship  were  not  adequate  to  examine  this 
possibility.  In  any  case,  it  seems  obvious  that  these  features  are  a  natural  phenomenon 
associated  with  low-grazing  angle  backscatter  from  the  sea  surface. 
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Figure  1 1 .  Space/time  images  of  ao(HH)  and  ao(W)  along  with  their  associated 
velocities.  Note  the  dramatically  larger  cross  section  and  velocities  for  FIH  polarization 
inside  the  pronounced  features.  The  bottom  panel  shows  ship  and  wind  parameters. 

To  document  this  farther,  and  to  assure  that  these  results  are  not  a  product  of  faulty 
radar  calibration,  we  show  time  series  of  ao(FlH)  and  ao(W)  at  an  incidence  angle  of 
88.75”  for  two  different  time  periods  near  the  Philippine  Islands  in  Figure  12.  This 
incidence  angle  was  chosen  because  it  corresponds  to  the  center  of  the  beam  from  the 
two-foot-diameter  parabolic  antennas  that  were  used  to  collect  these  data.  This  angle 
therefore  corresponds  to  the  largest  signal-to-noise  ratio.  One  of  the  time  periods  shown 
in  Figure  12  corresponds  to  the  image  of  Figure  1 1  while  the  other  was  obtained  well 


away  from  any  such  anomalous  backscattering  features.  Clearly  ao(HH)  is  smaller  than 
Oo(VV)  when  the  ocean  surface  is  disturbed  only  by  the  wind  but  becomes  larger  than 
ao(VV)  in  the  anomalous  region. 
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Figure  12.  Documentation  of  variable  polarization  ratios  at  an  incidence  angle  of  88.75“ 
+-  0.25°.  Top  row:  Solid  line  is  10  times  the  wind  speed,  circles  are  azimuth  angle  with 
respect  to  the  direction  from  which  the  wind  eomes,  X-  Bottom  row:  Measured  ao(HH) 
(red  asterisks)  and  ao(W)  (black  circles)  at  the  times  of  the  plots  in  the  top  row. 

Clearly  Bragg  seattering  cannot  account  for  these  features.  In  the  next  section,  we 
explore  some  ideas  that  have  been  put  forward  to  explain  these  non-Bragg  features  and 
take  the  first  steps  toward  a  method  for  including  the  most  important  ones  in  the  model. 

Other  Scattering  Mechanisms 


The  faet  that  the  multiseale  model  does  not  prediet  the  large  values  of  ao(HH)  that  are 
observed  at  incidence  angles  above  45°  in  all  look  direction  suggests  that  scattering 
processes  other  than  the  Bragg  scattering  processes  described  by  the  standard  composite 
surface  theory  may  be  occurring  at  high  incidence  angles.  This  has  been  noted  many 
times  in  the  past  by  a  variety  of  authors  (See  Churyumov  et  ah,  2002,  for  instance).  In 
fact,  the  data  suggest  that  more  than  one  of  these  non-standard  scattering  mechanisms 
may  be  at  work.  As  discussed  further  below,  many  of  the  non-Bragg  processes  that  have 
been  proposed  over  the  years  involve  breaking  waves;  Bragg  scattering  from  bound, 
tilted  waves  proposed  by  Plant  (1997)  also  require  breaking,  or  at  least  crumpling,  waves. 
However,  the  fact  that  ao(HH)  is  higher  than  predicted  in  the  downwind  look  direction 
(See  Figures  6-8)  is  difficult  to  explain  using  only  a  mechanism  tied  to  breaking  waves. 


Such  wave-breaking  regions  exist  primarily  near  wave  erests  or  on  the  front  faee  of 
waves  and  their  effeets  should  be  mueh  weaker  when  the  antenna  is  looking  downwind, 
espeeially  at  very  high  ineidenee  angles. 

Many  possibilities  for  the  additional  backseatter  at  HH  polarization  have  been 
proposed  over  the  years.  These  include  bound,  tilted  waves  due  to  the  breaking  of  short 
gravity  waves,  speeular  baekseatter,  wedges,  spray,  and  multipath  seattering  from  steep 
features. 

Plant  (1997)  proposed  that  some  features  of  HH  backseatter  at  ineidenee  angles 
greater  than  45“  ean  be  explained  by  Bragg  seattering  but  from  the  rough  patehes 
produeed  by  breaking  short  gravity  waves  that  travel  with  the  the  gravity  waves  rather 
than  from  freely  propagating  eentimetric  surface  waves.  In  partieular,  observations  of 
larger  Doppler  shifts  in  HH  than  in  VV  backseatter  and  of  low  eoherence  between  HH 
Doppler  shifts  and  dominant  oeean  wave  orbital  veloeities  when  looking  into  the  wind  at 
large  ineidenee  angles  ean  be  explained.  The  patehes  are  assumed  to  be  strongly  tilted 
toward  the  antenna  under  these  conditions  and  this  tilt  increases  HH  Bragg  scattering 
more  than  VV.  The  existenee  of  bound,  tilted  waves  in  wind-generated  and  breaking 
wave  systems  has  been  well  substantiated  by  wave  tank  experiments  (Plant  et  ah,  1999a, 
1999b,  2004;  Rozenberg  et  ah,  1999;  Lee  et  ah,  1999).  The  skewness  of  sea  surfaee 
slope  probability  distributions  ean  also  be  explained  on  the  basis  of  such  waves  (Plant, 
2003a).  Nevertheless,  bound,  tilted  waves  do  not  seem  to  aeeount  for  all  of  the  exeess 
baekseatter  (over  standard  Bragg/eomposite)  for  HH  polarization  at  high  ineidenee 
angles.  Bragg  seattering  from  bound,  tilted  waves  eannot  aeeount  for  large 
ao(HH)/ao(VV)  ratios  exeept  momentarily  (Plant,  2003b).  It  also  eannot  aeeount  for 
ao(HH)  being  too  large  in  the  downwind  look  direction.  Furthermore,  a  limit  exists  on 
how  strong  the  seattering  from  bound,  tilted  waves  ean  be  sinee  their  effeets  are  not 
observed  in  VV  baekseatter.  This  limit  is  suffieiently  low  to  ensure  that  bound,  tilted 
waves  cannot  explain  all  of  the  exeess  backseatter  even  in  the  upwind  direction. 

Specular  backseatter  has  frequently  been  invoked  as  a  means  for  inereasing  the 
ao(HH)/ao(VV)  polarization  ratio  at  large  ineidenee  angles  (Chubb,  1999;  Kudryatzev, 
2005).  While  the  assumption  that  suffieient  speeular  points  exist  nearly  perpendicular  to 
the  high-incidenee-angle  incident  field  to  produee  signifieant  backseatter  is  arguable,  this 
meehanism  clearly  eannot  produce  polarization  ratios  mueh  larger  than  one,  as  observed 
above. 

Wedge  seattering,  whieh  is  seattering  from  sharp  erests  of  signifieant  extent  in  the 
horizontal,  has  been  suggested  as  a  means  of  inereasing  high-incidence  angle  baekseatter 
at  HH  polarization  (Lyzenga  et  ah,  1983).  Previous  work  on  baekseatter  from  these 
struetures  had  been  eontradietory.  Kalmykov  and  Pustovoytenko  (1976)  suggested  that  it 
produeed  a  polarization  ratio  greater  than  one  at  high  ineidenee  angles  while  Lewis  and 
Olin  (1980)  said  it  produeed  polarization  ratios  less  than  one.  Lyzenga  et  al.  showed  that 
the  polarization  ratio  was  above  one  at  ineidenee  angles  higher  than  a  level  set  by  the 
wedge  angle  and  below  one  for  lower  ineidenee  angles.  Thus,  while  these  struetures  may 
eontribute  somewhat  to  the  baekseatter,  they  eannot  explain  our  observations  of  ao(HH)/ 
ao(VV)  >  1  at  very  high  ineidenee  angles. 

One  exeeption  to  this  statement  is  for  ineident  radiation  that  is  strikes  the  wedge  from 
above,  that  is,  perpendieular  to  its  base.  This  eonfiguration  is  most  likely  to  oeeur  for 
high  ineidenee  angles  when  the  “wedge”  is  the  tip  of  a  plunging  breaking  wave. 


Churyumov  et  al.  (2002)  showed  that  under  these  eonditions  the  baekseatter  may  have  a 
polarization  ratio  above  one.  Three-dimensional  numerieal  ealeulations  by  Li  and  West, 
(2006)  indieated  the  same  phenomenon.  The  density  of  plunging  breakers  on  the  open 
oeean  is  likely  to  be  rather  low  under  normal  conditions,  however,  so  this  mechanism 
may  contribute  relatively  little  to  the  mean  baekseatter  in  most  cases. 

Spray  above  the  water  surface  that  is  produced  by  actively  breaking  waves  is  another 
possible  source  of  baekseatter.  Bounces  from  the  surface  of  the  water  to  the  spray 
droplets  will  serve  to  augment  the  direct  baekseatter  from  the  drops.  Both  this  and 
possible  flattening  of  the  drops  could  produce  HH  baekseatter  that  is  larger  than  VV 
baekseatter.  However,  Plant  (2003b)  showed  that  if  spray  is  to  account  for  the  increased 
baekseatter  at  HH  polarization  and  high  incidence  angles,  it  must  increase  with  wind 
speed  much  more  slowly  than  is  generally  observed.  On  the  other  hand.  Plant  et  al. 

(2006)  later  showed  that  while  large  Doppler  shifts  due  to  baekseatter  from  spray  can  be 
observed  in  a  wind  wave  tank,  the  spray  responsible  for  the  baekseatter  is  located  too 
close  to  the  surface  to  be  measured  with  in  situ  instruments.  Thus  the  question  of  the  role 
of  spray  in  high-incidence  baekseatter  is  still  open. 

Multipath  scattering  from  steep  surface  features  is  perhaps  the  most  likely  source  of 
this  observation.  Baekseatter  from  breaking  waves  has  been  studied  for  many  years  and 
clearly  shows  that  ao(HH)  can  exceed  ao(VV)  during  breaking  (Lee  et  al,  1998;  Fuchs  et 
al.,  1999;  Sletten  et  al.,  2003) .  Most  researchers  have  invoked  some  type  of  multiple 
scattering  phenomena  to  explain  this  behavior  (Wetzel,  1986;  Trizna,  1997;  Holliday  et 
al.,  1998;  West,  1999,  2002;  Lee  et  al,  1997,  1998).  Wetzel  (1986)  seems  to  have  been 
the  first  to  suggest  the  importance  of  this  type  of  scattering.  He  proposed  that  the 
baekseatter  resulted  from  a  bounce  of  the  incoming  ray  from  the  water  surface  in  front  of 
a  bore-like  front  resulting  from  breaking.  He  modeled  this  as  a  quarter  cylinder  rising  out 
of  a  fiat,  tilted  sea  surface  and  multiplied  the  standard  cylinder  cross  section  by  a 
“proximity  factor”  to  account  for  the  presence  of  the  flat  surface.  His  results,  while 
including  Brewster  damping  of  vertically  polarized  rays  at  the  first  bounce,  indicate  that 
VV  baekseatter  will  be  stronger  than  HH  baekseatter  for  large  incidence  angles. 

Lee  et  al.  (1997)  and  Trizna  (1997)  both  proposed  that  the  structures  producing 
baekseatter  during  breaking  could  be  better  modeled  by  a  dihedral  in  which  two  water 
surfaces  at  right  angles  produce  a  double-bounce  baekseatter.  Both  suggested  that 
Brewster  damping  of  VV  polarized  incident  radiation  at  the  initial  bounce  would  reduce 
its  baekseatter  level  compared  to  that  of  HH  polarization.  Lee  et  al.  carried  out 
experiments  that  verified  this.  Trizna  concentrated  on  the  interference  of  the  ray  directly 
incident  on  the  vertical  surface  and  that  reflected  from  the  horizontal  surface.  Since  both 
of  these  rays  cannot  be  reflected  to  the  antenna  simultaneously,  however,  it  is  unclear 
why  this  interference  is  a  concern. 

While  dihedrals  are  clearly  idealized  scattering  structures,  numerical  studies  show  that 
similar  multiple  scattering  phenomena  occur  for  a  breaking  wave  at  and  immediately 
after  breaking  (West,1999,  2002;  Fuchs  et  ah,  1999).  These  studies  also  show  that  the 
exact  nature  of  this  multiple  scattering  depends  critically  on  the  shape  of  the  surface  (Lee 
et  al,  1998;  Fuchs  et  al.,  1999;  West,  2002;  Sletten  et  al.,2003).  Furthermore,  multiple 
scattering  also  occurs  for  strong  breaking  in  the  very  turbulent  region  produced  by 
breaking,  as  evidenced  by  the  cross  polarized  signal  and  ao(HH)  greater  than  ao(VV) 

(Lee  et  ah,  1998).  Therefore,  a  statistical  model  of  baekseatter  from  breaking  waves  in 


which  the  scattering  is  due  to  multiple  scattering  from  randomly  oriented  dihedrals  has 
been  proposed  (Lee  et  ah,  1998)  and  may  not  be  a  bad  approximation  to  reality. 

Water  Dihedrals 

As  a  first  step  in  this  proeedure,  we  ealeulate  the  baekseattering  eross  seetion  of 
dihedrals  formed  by  two  perpendieular  water  surfaees.  The  geometry  is  sketched  in 
Figure  13.  We  will  take  a  geometric  optics  approach  and  ignore  diffraction  except 
insofar  as  it  enters  into  the  determination  of  the  eross  seetion.  Thus,  we  eonsider  only 
ineident  and  baekseattered  rays  lying  in  the  plane  of  ineidenee  as  shown  in  Figure  12. 
We  will  let  the  eleetrie  field  veetor  lie  at  an  angle  (p  to  the  plane  of  ineidenee. 

Plane  of  Incidence 

/ 


Figure  13.  Diagram  of  the  dihedral  seattering  geometry.  Sides  of  the  dihedral  are  ay  and 
ah,  0g  is  the  grazing  angle,  and  the  eleetrie  field  E  makes  an  angle  (p  with  the  plane  of 
incidence. 

As  usual  for  simple  targets,  the  eross  seetion  a  is  given  by 

a  =  Ap|R|^D  (1) 

where  Ap  is  the  projeeted  area,  R  is  the  refieetivity,  and  D  is  the  direetivity.  We  will  let 
the  width  of  the  dihedral  by  b.  Then 

Ap=  b(aisin0g+ a2COS  0g),  ai  =  min(ah,aveot0g);  a2  =  min(av,ahtan0g) 

R=  [Rv(0g)exp(-i7i/2)eos(p  +  Rh(0g)sin(p][Rv(7i/2-0g)  exp(-i7i/2)eos(p  +  Rh(7i/2-0g)sin(p]  (2) 
D  =  4  TiA^k^ 

where  Rv  and  Rh  are  Fresnel  refieetion  eoeffieients  for  vertieal  and  horizontal 
polarization,  respeetively,  X  is  mierowave  length,  and  ay  and  ah  are  the  lengths  of  the 
vertieal  and  horizontal  planes  of  the  dihedral.  The  faetor  exp(-i7i/2)  in  the  refieetivity 
aeeounts  for  the  faet  that  Ry  is  a  ratio  of  magnetie  fields  while  Rh  is  a  ratio  of  eleetrie 
fields.  For  perfeet  refieetivity,  Ry  =  1  and  Rh  =  -1  so  the  maximum  eross  seetion  at  0g  = 
45°  with  ah  =  a2  =  a  is 

Omax  =  87ia^b^/k^ 


(3) 


as  required.  The  eross  seetion  for  a  water  dihedral  at  X-band  (^  =  3  em)  is  plotted  versus 
grazing  angle  in  Figure  14  for  b  =  1  em  and  various  values  of  a.  To  produee  these  eurves, 
Eq.(l)  has  been  eonvolved  with  a  spreading  funetion  eorresponding  to  the  vertieal 
beamwidth  of  the  vertical  plate; 

f(0g)  =  exp  {-[4  ln2  0g^  }/N  (4) 

where  N  is  a  normalizing  factor  to  make  the  integral  of  f  over  0g  be  one. 
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Figure  14.  Cross  sections  of  dihedral  comer  reflectors  of  various  sizes  as  a  function 
grazing  angle.  The  curves  with  diamonds  are  HHA^V  polarization  ratios  in  dB.  The 
width  of  the  dihedral  is  1  cm.  ai  is  the  length  of  the  horizontal  plane  while  a2  is  the 
length  of  the  vertical  plane. 

Several  features  of  these  plots  are  of  interest.  First,  while  the  cross  sections  obviously 
increase  with  dihedral  size  at  45°,  this  is  not  tme  at  the  grazing  angles  of  more  interest  to 
sea  clutter.  In  particular,  for  grazing  or  incidence  angles  near  5°,  cross  sections  do  not 
monotonically  increase  with  dihedral  size.  Nevertheless,  these  levels  alone  do  not 
determine  the  magnitudes  of  backscattering  cross  sections  since  the  temporal 
and  spatial  densities  of  the  dihedrals  will  also  affect  the  values.  Second,  peaks  of  the 
cross  sections  at  0°  and  90°  grazing  angles  are  due  to  specular  reflections  from  the 
vertical  and  horizontal  planes  of  the  dihedrals,  respectively.  They  will  clearly  affect 
cross  sections  not  only  near  grazing  but  also  near  nadir.  The  nadir  cross  sections  may 
limit  the  densities  of  these  types  of  scatterers  that  we  can  choose  since  backscattering  at 
nadir  is  well  understood.  Third,  note  that  the  HHAX^  polarization  ratio  is  greater  than 
one  only  at  grazing  angles  above  a  few  degrees.  Since  this  is  often  higher  than 
the  nominal  grazing  angles  used,  for  instance,  from  ships,  observations  of  the  polarization 
ratios  will  be  affected  by  the  tilting  of  the  dihedrals  by  the  long  waves  on  which  they  ride. 
Finally,  note  that  for  grazing  angles  between  about  5°  and  20°,  the  polarization  ratio 
increases  with  dihedral  size.  This  effect  will  not  be  changed  by  summing  over  many  such 
features.  Thus  we  suspect  that  larger  observed  polarization  ratios  will  correspond  to 
larger  dihedral  features  on  the  surface  and  therefore  to  more  energetic  breaking  events. 


Conclusion 


Measurements  of  the  mean  normalized  radar  eross  seetion  of  the  sea  taken  over  many 
years,  from  many  different  platforms,  and  at  many  different  ineidenee  angles  have  shed 
light  on  possible  baekseattering  meehanisms  of  the  oeean  surfaee.  The  measurements 
show  that  ao(VV)  for  a  sea  surfaee  disturbed  only  by  wind  is  generally  well  explained  by 
the  multiseale  model  (Plant,  2002)  over  the  entire  range  of  ineidenee  angles  from  0°  to 
89°.  This  model  is  simply  a  slightly  modified  version  of  the  standard  eomposite  surfaee 
theory  of  Wright  (1968)  and  Bass  et  al.  (1968).  ao(HH),  on  the  other  hand,  begins  to  be 
higher  than  predieted  by  this  model  at  ineidenee  angles  above  about  45°  but  remains 
below  ao(W)  under  normal  eonditions.  While  at  very  high  ineidenee  angles,  both  eross 
seetions  inerease  with  wind  speed  as  in  standard  seatterometry  and  ao(VV)  displays  the 
standard  seeond-harmonie  dependenee  on  azimuth  angle,  ao(HH)  shows  a  very  different 
azimuth  angle  dependenee.  It  is  highest  looking  into  the  wind  and  lowest  looking 
downwind.  In  both  oases  the  magnitude  of  ao(HH)  is  larger  than  expeoted  from  the 
multiseale  model. 

Observations  of  anomalous  baoksoatter  from  the  oeean  that  depart  from  this  model 
have  been  obtained  in  the  Philippine  Sea.  While  the  oause  of  this  anomalous  baoksoatter 
is  unknown,  it  olearly  shows  that  ao(HH)  oan  at  times  exoeed  ao(VV)  over  large  time  and 
spaoe  soales  for  low-grazing-angle  baoksoatter  from  the  oeean.  Taken  together,  these 
observations  show  that  other  soatterers  besides  short  wind-driven  waves  exist  on  the 
oeean  surfaee.  Several  different  types  of  these  soatterers  probably  exist.  They  inolude 
bound,  tilted  waves  from  “orumpling”  short  gravity  waves,  jets  from  plunging  breaking 
waves,  spray,  and  multipath  surfaoes  most  oonveniently  desoribed  as  water  dihedrals. 
These  types  of  soatterers  begin  to  manifest  themselves  when  their  normalized  radar  eross 
seetion  exoeeds  the  standard  Bragg  soattering  from  wind  waves.  This  most  oommonly 
ooours  for  HH  polarization  at  high  ineidenee  angles  but  may  also  be  oaused  or 
exaggerated  by  anomalous  eonditions  suoh  as  shear  ourrents  as  we  will  show  in  Part  2. 

Sinoe  multipath  surfaoes  are  likely  to  be  the  dominant  types  of  these  other  soatterers 
at  very  high  ineidenee  angles,  we  derived  a  geometrioal  optios  expression  for  their  eross 
seetion  versus  grazing  angle.  This  showed  that  suoh  dihedrals  must  be  tilted  up  toward 
the  inooming  signal  in  order  to  produoe  ao(HH)  /ao(VV)  polarization  ratios  that  exoeed 
one.  Still  to  be  determined  are  the  distribution  of  orientations  of  the  dihedrals  and  their 
density  on  the  oeean  surfaee  under  different  oiroumstanoes.  A  limitation  on  these 
variables  will  be  the  known  density  of  speoular  soatterers  at  nadir. 
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We  report  measurements  of  the  surface  signatures  of  internal  waves  with  dual- 
polarized,  coherent,  X-band  radars  mounted  on  three  different  ships  and  an  airplane. 

Two  sets  of  shipboard  measurements  were  made  in  the  South  China  Sea  while  the 
airborne  data  and  the  third  set  of  shipboard  data  were  collected  in  the  Atlantic  Ocean  off 
the  New  Jersey  coast.  In  the  shipboard  measurements  in  the  South  China  Sea,  internal 
waves  generally  increased  the  backscattering  cross  section  near  the  peaks  of  the  internal 
waves  with  no  detectable  decrease  afterward.  A  very  few  times,  however,  the  signatures 
did  show  increases  followed  by  decreases  in  cross  section.  The  peak  of  the  cross  section 
signature  shifted  its  location  relative  to  the  internal  wave  crest  depending  on  the 
maximum  strain  rate  of  the  internal  wave.  In  the  shipboard  measurements  off  the  New 
Jersey  coast,  the  range  resolution  of  the  radar  was  sufficiently  small  to  allow  us  to  detect 
significant  modulations  of  surface  waves  on  the  order  of  15  to  30  meters  long  by  the 
internal  waves.  On  one  low-wind-speed  day,  simultaneous  and  collocated  airborne  and 
shipboard  measurements  were  obtained  off  the  New  Jersey  coast.  Significant  differences 
were  seen  between  internal  wave  signatures  at  the  very  high  incidence  angles  of  the 
shipboard  measurements  and  those  at  the  moderate  incidence  angles  of  the  airborne  data. 
At  the  high  incidence  angles,  the  cross  section  for  horizontally  polarized  transmit  and 
receive  signals,  ao(HH),  regularly  exceeded  that  for  vertically  polarized  transmit  and 
receive  signals,  ao(VV),  near  the  internal  wave  crest  by  5  to  15  dB.  At  the  more 
moderate  incidence  angles  observed  from  the  aircraft,  maximum  values  of  ao(HH)  and 
ao(VV)  are  more  nearly  equal,  with  ao(HH)  often  being  comparable  to  ao(VV)  and  only 
occasionally  exceeding  it.  The  intensity  of  modulation  of  the  cross  section  caused  by 
internal  waves  observed  from  the  plane  depended  little  on  the  direction  of  observation. 

In  all  cases,  internal  wave  surface  signatures  became  less  visible  with  increasing  wind 
speed,  being  very  difficult  to  observe  at  9  m/s. 

Introduction 

When  internal  waves  propagate  on  density  interfaces  below  the  ocean  surface,  they 
produce  surface  currents  that  modulate  surface  waves,  thus  producing  detectable 
signatures  in  microwave  backscatter  from  the  surface.  These  signatures  are  easily 
observed  by  synthetic  aperture  radar  and  a  very  large  number  of  images  of  internal  wave 
surface  signatures  have  been  obtained  from  orbiting  spacecraft  (Apel,  2004,  Liu  et  ah, 
2008).  Recently,  sea  return  to  shipboard  radars  has  begun  to  be  recorded  and  internal 
wave  surface  signatures  in  this  return  have  been  examined  (Liu,  et  ah,  2008,  Chang  et  ah, 
2008). 

Studies  aimed  at  understanding  the  mechanisms  by  which  internal  waves  modulate 
surface  waves  to  produce  microwave  signatures  have  been  carried  out  over  many  years. 


Among  the  experimental  studies  are  the  seminal  observations  of  Apel  et  al.  (1975)  from 
the  Earth  Resourees  Teehnology  satellite,  the  1972  work  of  Hughes  and  Grant  (1978),  the 
1983  Joint  Canada-US  Oeean  Wave  Investigation  Project  (JOWIP,  also  called  the 
Georgia  Straits  experiment)  (Hughes  and  Dawson,  1988),  the  1984  SAR  Internal  Wave 
Signature  Experiment  (SARSEX)  (Hughes  and  Gasparovic,  1988),  the  1987  Eock  Einnhe 
experiment  (Hogan  et  ah,  1996),  and  the  1995  Coastal  Ocean  Probing  Experiment 
(COPE)  (Kropfli  et  ah,  1999).  Theoretical  studies  into  the  generation  of  internal  wave 
microwave  surface  signatures  include  Alpers  (1985),  Eyzenga  and  Bennett  (1988), 
Thompson,  1988,  and  Kudryavtsev  et  al.  (2005). 

Nevertheless,  significant  gaps  still  exist  in  our  understanding  of  how  the  surface 
currents  associated  with  internal  waves  generate  the  signatures  observed  in  microwave 
backscatter.  While  the  role  of  surface  current  gradients  in  straining  surface  waves  is 
clear,  it  is  still  not  known  how  often  these  current  gradients  simply  modulate  surface 
waves  and  how  often  they  cause  surface  waves  to  break  thus  generating  the  short 
microwave  scattering  elements  on  the  surface.  Because  both  transmitting  and  receiving  a 
horizontally  polarized  microwave  signal  (HH  backscatter)  causes  sea  return  to  be  more 
sensitive  to  breaking  surface  waves  than  transmitting  and  receiving  vertically  polarized 
signals  (VV  backscatter),  comparison  of  internal  wave  signatures  at  the  two  polarizations 
offers  the  possibility  of  investigating  the  prevalence  of  breaking  waves  in  the  signatures 
(Plant,  1997;  Eiu  et  ah,  1998;  Hwang  et  ah,  2008).  Unfortunately,  very  little  work  has 
been  done  in  which  HH  and  VV  polarized  signals  have  been  used  at  nearly  the  same  time 
to  observe  internal  wave  signatures,  especially  at  moderate  incidence  angles  (Kropfli, 
1999;  Churyumov,  2002).  Even  less  work  has  been  done  to  produce  well-calibrated 
NRCS  values  near  internal  waves.  Similarly,  the  determination  of  the  dependence  of  the 
microwave  signatures  on  azimuth  angle  has  been  little  studied.  Einally,  in  most 
situations,  it  is  difficult  to  determine  precisely  the  location  of  the  surface  signature  with 
respect  to  an  internal  wave  crest. 

In  this  paper,  we  address  these  issues.  We  report  the  effect  of  internal  waves  (IWs) 
on  the  normalized  radar  cross  section  of  the  sea,  NRCS  or  Oq,  and  on  the  surface  velocity, 
which  for  microwaves  is  really  the  velocity  of  the  scatterers  of  the  microwave  radiation. 
All  measurements  reported  here  were  made  with  HH  and  VV  returns  measured  within 
seconds  of  each  other  (although  in  one  experiment  a  switch  failed  so  we  only  obtained 
VV  signals)  and  all  NRCS  values  are  calibrated.  While  the  shipboard  measurements 
were  generally  made  with  the  radar  antennas  looking  either  in  the  direction  of  internal 
wave  propagation  or  against  it,  the  aircraft  measurements  were  made  at  many  different 
azimuth  angles  so  the  dependence  of  signature  intensities  on  antenna  look  direction  could 
be  studied.  Einally,  by  comparing  the  shipboard  radar  measurements  with  currents 
measured  simultaneously  by  acoustic  Doppler  current  profilers  (ADCPs)  on  the  same 
ships,  the  location  of  the  surface  signature  relative  to  IW  crests  could  be  determined. 

Because  our  measurements  were  made  both  in  the  South  China  Sea  and  in  the 
Atlantic  Ocean  off  the  New  Jersey  coast,  two  very  different  types  of  internal  waves  were 
studied.  In  the  South  China  Sea,  surface  signatures  of  very  non-linear  internal  solitons, 
which  traveled  in  widely  spaced  trains,  were  observed.  We  also  found  that  in  the  South 
China  Sea,  regular  nearly  sinusoidal  trains  of  smaller  amplitude  internal  waves  were  very 
frequently  observed.  These  smaller  waves  were  very  similar  to  the  internal  waves 


encountered  off  the  New  Jersey  coast.  We  first  discuss  signatures  of  the  non-linear 
solitons  then  the  smaller  sinusoidal  wave  trains. 


Shipboard  Measurements  in  the  South  China  Sea 

Two  sets  of  shipboard  measurements  were  made  in  the  South  China  Sea.  In  2005,  our 
X-band,  coherent  radar  called  RiverRad  was  mounted  on  the  RW  Revelle  and  took 
measurements  in  the  South  China  Sea  from  April  18  to  May  14.  These  measurements 
were  all  taken  in  deep  water.  In  2007,  we  operated  RiverRad  on  two  cruises,  one  from 
April  24  to  May  13  was  in  the  relatively  shallow  water  of  the  western  shelf  of  the  South 
China  Sea  while  the  second  from  May  13  to  May  20  was  in  deep  water.  On  all  these 
cruises,  the  parabolic,  pencil-beam  antennas  were  pointed  approximately  in  the  direction 
of  the  ship’s  heading.  Figure  1  shows  RiverRad  mounted  on  these  two  ships  while  Figure 
2  indicates  where  in  the  South  China  Sea  significant  non-linear  solitons  were 
encountered. 


Figure  1.  a)  The  parabolic  antennas  of  RiverRad  are  seen  mounted  on  the  ICV  Revelle  in 
2005.  b)  RiverRad  mounted  on  the  Taiwanese  ship  ICV  Ocean  Researcher  1  in  2007. 


Symbols:  x  =  2005,  deep;  o  =  2007,  shelf;  +  =  2007.  deep 


Figure  2.  Locations  at  which  non-linear  internal  solitons  were  encountered  in  the  South 
China  Sea.  Shallow  water,  2007  -  circles;  deep  water,  2007  -  +;  deep  water,  2005  -  X. 


Figure  3  shows  the  types  of  data  collected  on  these  cruises.  The  standard  shipboard 
measurements  of  wind  speed  and  direction,  ship  speed  and  direction,  and  currents  from 
ADCPs  were  collected  and  recorded.  From  the  latter  measurements,  currents  were 
extracted  from  the  nearest  good  bin  to  the  surface  and  resolved  into  components  in  and 
perpendicular  to  the  ship  heading  as  shown  in  the  third  panel  from  the  top  of  Figure  3.  In 
2005,  the  nearest  good  bin  to  the  surface  was  17  m  below  the  surface;  in  2007,  it  was  52 
m  deep.  In  addition,  latitude,  longitude,  time,  pitch,  and  roll  were  recorded  during  the 
cruise.  The  latter  two  variables  were  used  in  correcting  the  incidence  angle  of  the  radar 
for  ship  motion.  We  also  had  a  GPS  receiver  mounted  on  the  RiverRad  antennas  and  its 
speed  readings  were  used  to  correct  the  Doppler  shifts  observed  by  RiverRad.  In  addition 
to  the  ship  speed  from  the  GPS  shown  in  Figure  3,  the  speed  indicated  by  a  pitot  tube  on 
the  ship’s  hull  was  recorded.  The  difference  between  the  two  measured  speeds  was 
another  measure  of  current  in  the  heading  direction. 
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Figure  3.  An  example  of  data  collected  from  the  R/V  Revelle  in  2005  averaged  over  50 
sec.  The  top  panel  shows  radar  surface  (scatterer)  velocities  in  the  direction  of  the  ship’s 
heading  as  functions  of  distance  in  the  heading  direction  and  time.  Distances  are 
wrapped  at  4  km.  The  curve  of  black  dots  is  the  ship’s  location.  The  second  panel  down 
shows  a  similar  plot  of  normalized  radar  cross  section  in  the  direction  of  the  ship’s 
heading.  The  microwave  measurements  were  made  only  at  W  polarization  in  2005. 

The  third  panel  down  shows  the  ship  speed  in  green,  the  wind  speed  in  red,  four  times  the 
ADCP  current  at  17  m  depth  along  (black)  and  perpendicular  (blue)  to  the  ship’s  heading. 
Currents  in  the  heading  direction  are  positive  if  they  come  toward  the  ship  and  currents 
perpendicular  to  the  heading  are  positive  to  starboard.  Finally  the  bottom  panel  shows 


compass  readings  of  the  ship’s  heading  (black),  ship  track  (green),  and  wind  direetion 
(from,  red). 

Mierowave  data  consisted  of  the  standard  output  of  RiverRad,  whieh  has  been 
deseribed  previously  elsewhere  (Plant,  2005b).  Briefly,  RiverRad  is  a  coherent,  9.36 
GHz  radar  that  reeeives  signals  from  a  series  of  range  bins  sufficiently  rapidly  to 
calculate  a  Doppler  spectrum  at  every  bin.  Averages  of  these  speetra  over  a  complete 
record  were  stored  for  each  range  bin.  In  2005  a  reeord  was  50  seeonds  long,  in  2007  it 
was  30  minutes  long.  In  addition,  the  first  three  moments  of  the  Doppler  spectrum,  mean 
power,  Doppler  offset,  and  Doppler  bandwidth  are  stored  for  every  scan,  each  of  which 
was  a  approximately  365  msec  long.  Alternate  seans  were  eollected  with  HH  and  VV 
polarizations.  Therefore  a  scan  at  each  polarization  was  365  msec  long  and  scans 
occurred  every  730  msee.  The  range  resolution  was  15  m. 

From  the  three  reeorded  moments  of  the  speetrum,  we  obtained  the  normalized  radar 
eross  seetion  (NRCS  or  Oo,  a  measure  of  surfaee  roughness),  the  mean  surface  veloeity, 
and  the  spread  in  scatterer  veloeities  for  each  range  bin  averaged  over  a  scan  time  (365 
msee).  Images  of  the  surface  veloeity  are  shown  in  the  top  panel  of  Figure  3  while 
images  of  Oo  are  shown  in  the  seeond  panel  from  the  top.  These  images  are  space-time 
plots  of  the  return  to  the  radar  along  the  antenna-look  direction,  the  direction  of  the  ship’s 
heading.  They  were  smoothed  over  50  see.  The  swaths  wrap  from  top  to  bottom  at  4  km 
so  the  vertieal  axis  is  distance  along  the  ship’s  heading  from  the  start  of  the  run,  modulo 
4.  The  black  dotted  eurves  in  the  images  show  the  loeation  of  the  ship.  The  radar  return 
starts  about  250  m  in  front  of  the  ship  and  goes  out  to  almost  2  km.  The  antennas  were 
18  m  above  the  water  surfaee  so  the  incidence  angle  ranged  from  85.9'’  to  89.5°.  Note 
that  this  type  of  display  plaees  the  radar  return  at  the  loeation  it  has  for  a  stationary 
observer  if  the  ship’s  heading  is  eonstant. 

Surface  signatures  of  two  types  of  internal  waves,  small  nearly  sinusoidal  waves  and  a 
larger  soliton  wave,  are  clearly  present  in  both  the  velocity  and  cross  section  measured  by 
the  radar.  As  indicated  by  the  sign  of  the  ADCP-measured  internal  soliton  current  along 
the  heading  and  the  lack  of  internal  soliton  current  perpendicular  to  the  heading,  the  ship 
was  travelling  direetly  opposite  the  direetion  of  propagation  of  the  large  internal  wave. 
The  train  of  smaller  sinusoidal  internal  waves  with  a  wavelength  of  less  than  800  m  is 
observed  in  both  the  surfaee  eurrents  and  eross  sections  in  front  of  the  large  soliton.  Both 
the  surfaee  veloeity  and  the  cross  section  of  the  smaller  waves  inerease  significantly  over 
the  large  internal  wave.  Apparently  the  interaction  of  the  soliton  and  the  smaller  wave 
train  increases  the  amplitude  of  the  small  wave  train.  In  fact  the  small  wave  train  is 
difficult  to  observe  after  the  ship  has  passed  the  soliton,  indicating  that  it  may  have  been 
blocked  by  the  soliton. 

The  slope  of  these  surfaee  features  is  the  speed  of  the  internal  wave  if  the 
perpendicular  ADCP  current  is  zero.  If  the  internal  waves  do  not  propagate  directly 
toward  or  away  from  the  radar,  then  the  slope  of  the  feature  is  the  maximum  possible 
speed  of  the  internal  wave.  The  train  of  small  internal  waves  is  clearly  visible  in  the 
radar  return  but  is  barely  visible  in  the  currents  measured  at  17  m  depth.  This  makes  it 
difficult  to  determine  their  direction  of  propagation  except  that  a  component  of  their 
propagation  vector  comes  toward  the  ship  beeause  the  slope  is  negative.  We  ean  say  with 


certainty  that  the  train  of  small  internal  waves  is  moving  mueh  more  slowly  than  the 
soliton.  Furthermore,  their  speed  is  increasing  as  they  approach  the  erest  of  the  soliton. 

In  order  to  obtain  the  precise  loeation  of  the  inereases  in  eross  seetion  and  surface 
velocity  caused  by  the  soliton  with  respect  to  the  location  of  the  soliton,  their  effeets  must 
be  extrapolated  to  the  position  of  the  ship.  Carrying  out  this  extrapolation  in  the  images 
of  Figure  3,  clearly  shows  that  the  surfaee  signature  of  the  internal  wave,  in  both  cross 
section  and  velocity  is  nearly  directly  above  the  soliton  crest.  As  we  will  show  later,  we 
generally  found  that  in  deep  water,  the  surface  signature  maximized  near  the  crest  of  the 
internal  wave. 

Figure  4  shows  that  the  signature  does  not  always  maximize  direetly  above  the  soliton. 
The  data  shown  in  this  figure  are  from  the  2007  eruise  on  the  ORl  in  shallow  water.  The 
format  of  this  figure  is  the  same  as  that  of  Figure  3.  Again,  the  ship  travels  very  nearly 
opposite  the  direetion  of  soliton  movement.  Now,  however,  maxima  in  both  the  cross 
section  and  surface  veloeity  signatures  of  the  wave  occur  well  before  the  erest  of  the 
soliton  reaches  the  ship. 
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Figure  4.  Data  from  the  2007  eruise  of  the  ORl  on  the  shelf  of  the  South  China  Sea.  The 
format  is  the  same  as  that  of  Figure  3. 


Figure  5  summarizes  the  loeation  of  the  maximum  cross  section  with  respect  to 
soliton  wave  crests  for  all  of  the  data  we  eolleeted  in  2005  and  2007.  The  figure  plots  the 
distance  by  which  the  signature  maximum  leads  the  wave  erest  versus  the  maximum 
eurrent  produced  by  the  wave  at  17  or  53  m  depth,  both  normalized  by  the  width  of  the 
wave.  The  definition  of  width  was  taken  to  be  the  full  width  at  half  the  maximum 
eurrent.  The  figure  shows  that  the  signature  of  the  soliton  oceurs  farther  in  front  of  the 
erest  the  steeper  the  wave  is.  In  general  solitons  in  the  shallow  water  of  the  shelf  are 
steeper  than  in  deep  water  and  their  signature  leads  the  soliton  crest  by  greater  distances. 
In  only  one  ease,  the  x  enelosed  in  a  diamond,  did  this  behavior  not  oeeur. 


In  the  vast  majority  of  cases  that  we  examined,  the  response  of  the  cross  section  to  the 
surface  currents  set  up  by  the  internal  waves  was  an  increase  unaccompanied  by  any 
significant  decrease.  The  symbols  in  Figure  5  that  are  enclosed  in  squares  show  the  cases 
where  this  was  not  true  and  the  increase  was  followed  by  a  significant  decrease.  We  do 
not  believe  that  the  behavior  of  the  majority  of  our  data  is  due  to  the  fact  that  a  decrease 
in  cross  section  is  hidden  by  noise  as  suggested  by  da  Silva  et  al.  (1998)  for  SAR  imagery 
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Figure  5.  Distance  by  which  the  maximum  cross  section  leads  the  internal  wave  crest  as 
a  fraction  of  the  width  of  the  wave  versus  the  maximum  current  set  up  by  the  internal 
wave  divided  by  the  width  of  the  wave,  a  measure  of  the  maximum  strain  rate.  Symbols 
with  squares  or  a  diamond  around  them  are  discussed  in  the  text. 


Figure  6.  a)  Colored  curves  show  the  NRCS  measured  at  various  ranges  as  a  function  of 
time.  NRCS  values  at  different  ranges  are  offset  as  indicated  in  the  figure.  The  black 
curve  is  30  times  the  ADCP  current  in  the  heading  direction,  b)  Black  symbols  show  the 
range  dependence  of  the  NRCS  measured  at  the  time  indicated  by  the  vertical  red  line  in 
a.  The  green  curve  shows  the  noise  level  of  the  radar,  c)  The  red  curve  is  ship  speed,  the 
blue  curve  is  ship  heading  divided  by  100,  the  black  curve  is  four  times  the  ADCP 
current  in  the  heading  direction,  and  the  green  curve  is  the  ADCP  current  perpendicular 
to  the  heading. 
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Figure  7.  Same  as  Figure  6  but  for  April  28,  2005,  illustrating  a  ease  where  an  internal 
wave  caused  an  NRCS  increase  followed  by  a  decrease. 
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Figure  8.  Illustration  of  two  sequential  internal  wave  causing  different  responses  in  the 
NRCS.  The  top  panel  is  a  time  series  of  VV  NRCS  values  at  various  ranges  as  the 
internal  solitons  are  crossed.  The  second  panel  from  the  top  is  a  similar  plot  for  HH 
polarization.  In  both  these  panels,  NRCS  values  are  offset  at  various  ranges  as  shown  on 
the  figure.  The  third  panel  from  the  top  gives  the  polarization  ratios  at  various  ranges 
without  offsets.  The  bottom  panel  shows  the  ship  speed  in  red,  the  ship  heading  divided 
by  100  in  blue,  the  along-heading  ADCP  current  times  four  in  black,  and  the  cross¬ 
heading  ADCP  current  in  green. 


of  internal  waves.  Figure  6  illustrates  that  RiverRad’s  signals  are  generally  well  above 
the  noise  level.  Two  of  the  cases  where  increases  were  followed  by  decreases  are 
illustrated  in  Figures  7  and  8.  Figure  7  is  in  the  same  format  as  Figure  6  and  shows  an 
NRCS  increase  followed  by  a  decrease.  Note  the  presence  of  small  nearly  sinusoidal 


waves  before  the  internal  soliton  arrives  at  the  ship  and  their  absenee  afterward.  Figure 
8  illustrates  a  ease  from  2007  where  the  sinusoidal  waves  are  not  so  prominent.  Note  that 
the  inerease/deerease  oeeurs  for  both  HH  and  VV  polarizations  for  the  first  internal  wave 
but  does  not  oeeur  for  either  when  the  seeond  internal  wave  arrives.  For  both  internal 
waves,  the  polarization  ratio  exeeeds  0  dB  at  its  maximum  value. 

We  have  eompared  these  observations  with  the  output  of  the  standard  aetion  balanee 
model  of  short  wave  modulation  by  internal  wave  eurrents  (Thompson  and  Gasparovie, 
1986;  Thompson,  1988;  Lyzenga  and  Bennett,  1988;  Kudryatsev  et  ah,  2005).  Most 
researchers  agree  that  the  direct  modulation  of  centimeter-length  Bragg  scatterers  by  the 
internal  wave  surface  currents  is  very  small.  Therefore  the  consensus  is  that 
intermediate-scale  waves  are  modulated  by  the  internal  wave  currents  and  that  at  the 
location  where  these  maximize,  the  modulation  of  centimetric  wind  waves  or  the 
generation  of  small-scale  roughness  by  breaking  will  also  maximize.  Thus  for  the 
purposes  of  our  comparison,  we  looked  at  the  modulation  of  the  mean-square  slope  of 
waves  between  24  cm  and  6.67  m  long.  We  computed  using  the  action  balance 
techniques  outlined  by  Lyzenga  and  Bennett  (1988)  the  relative  modulation  of  these 
waves  by  internal  solitons  of  fixed  width  but  of  variable  amplitude.  We  did  this  for  three 
different  wind  directions  relative  to  the  internal  wave  propagation  direction,  0°,  90°,  and 
180°.  The  results  are  shown  in  Figure  9. 


IW  Current  Profile,  Constant  in  y  Dir;  Wind  Speed  =4  m/s,  Wind  Dir  =  180  deg  CW  from  +x  Dir 


IW  Current  Profile,  Constant  in  y  Dir;  Wind  Speed  =4  m/s,  Wind  Dir  =  90  deg  CW  from  +x  Dir 


IW  Current  Profile,  Constant  in  y  Dir;  Wind  Speed  =4  m/s,  Wind  Dir  =  0  deg  CW  from  +x  Dir 
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Figure  9.  a)  Surface  currents  for  internal  waves  of  six  different  amplitudes  with  the 
positions  of  the  maximum  mean  square  slope  of  short  gravity  waves  between  24  and  667 
cm  long  as  predicted  by  action  balance  equations  shown  as  circles,  b)  Mean-square  slope 
of  short  gravity  waves  relative  to  the  mean-square  slope  of  a  wave  field  undisturbed  by 
IW  currents  for  the  six  different  solitary  wave  currents  shown  in  a.  In  a  and  b,  the  4  m/s 
wind  blows  opposite  the  IW  propagation  direction,  c)  Same  as  a  but  wind  now  blows 
perpendicular  to  the  IW  propagation  direction,  d)  Same  as  b  but  wind  now  blows 
perpendicular  to  the  IW.  e)  Same  as  a  but  wind  now  blows  in  the  direction  of  the  IW 
propagation  direction,  f)  Same  as  b  but  wind  now  blows  in  the  direction  of  the  IW 
propagation  direction. 

We  found  that  when  the  wind  blew  against  the  internal  wave  (Figures  9a  and  b),  the 
maximum  modulation  moved  down  the  forward  face  of  the  internal  wave  as  its  amplitude 
increased  but  by  much  smaller  distances  than  we  observed.  When  the  wind  blew 
perpendicular  to  the  internal  wave  (Figures  9c  and  d),  the  maximum  modulation  moved 
much  more  in  line  with  our  observations.  When  the  wind  blew  in  the  direction  of  the 
internal  wave,  the  maximum  modulation  moved  down  the  forward  face  of  the  internal 
wave  somewhat  farther  than  the  majority  of  our  measurements  did  (Figures  9e  and  f). 
Figure  10  shows  histograms  of  wind  speeds  and  directions  relative  to  the  IW  propagation 
directions  for  the  times  of  our  measurements. 

For  all  three  of  the  wind  directions  that  we  modeled,  a  depression  of  the  relative 
modulation  occurred  behind  the  crest  of  the  internal  wave.  This  depression  was  generally 
small  for  low-amplitude  internal  waves  but  became  larger  as  the  internal  wave  steepness 
increased.  Our  data  do  not  support  the  idea  of  a  large  decrease  of  NRCS  being  associated 
with  large  internal  wave  steepnesses  (see  Figure  5).  Nevertheless,  the  small  decreases 
associated  with  smaller  amplitude  internal  solitons  may  explain  why  we  generally  did  not 
observe  decreased  NRCS  associated  with  the  surface  signatures.  If  the  decrease  was  as 
small  as  the  calculations  suggest,  we  may  not  have  been  able  to  detect  it. 
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Figure  10.  a)  Flistogram  of  wind  speeds  during  measurements  of  internal  waves,  b) 
Histogram  of  wind  directions  relative  to  the  IW  propagation  direction  during  the  same 
measurements. 


Airborne  and  Shipboard  Measurements  off  the  New  Jersey  Coast 

Measurements  in  the  Atlantic  Ocean  near  the  New  Jersey  shelf  break  were  made  both 
on  the  R/V  Endeavor  and  on  a  Cessna  Skymaster.  Figure  1 1  shows  RiverRad  mounted 
on  the  ship  (left)  and  a  similar  coherent  radar,  CORAR,  mounted  on  the  plane  (right). 

For  the  shipboard  measurements,  RiverRad’s  antennas  were  scanned  back  and  forth  in 
a  windshield-wiper  fashon  through  an  arc  80°  in  azimuth  centered  at  24°  higher  azimuth 
angle  than  the  heading.  The  antennas  were  mounted  35°  apart  in  azimuth;  one  was  W 
polarized  and  the  other  was  HH  polarized.  The  W  antenna  was  at  the  higher  azimuth 
angle.  The  antennas  swept  through  their  80°  arc  in  approximately  40  seconds.  The  range 
resolution  was  set  to  3.75  m. 

The  airborne  measurements  were  made  in  a  side-looking  mode.  The  W  and  HH 
antennas  were  mounted  one  above  the  other  looking  to  the  left  side  of  the  plane.  CORAR 
has  been  described  in  detail  by  Plant  et  ah,  2005a.  CORAR  is  very  similar  to  RiverRad 


Figure  11.  Feft;  RiverRad  mounted  the  R/V  Endeavor  in  2006.  Right:  CORAR  mounted 
on  the  Cessna  Skymaster  at  the  same  time. 


except  that  data  are  collected  and  stored  more  quickly  as  is  necessary  for  airborne 
operation.  A  scan  was  340  msec  in  length  and  the  first  three  spectral  moments  were 
saved  at  every  range  bin  on  each  scan  along  with  a  complete  spectrum  from  a  selected 
range  bin. 

Figure  12  shows  an  image  of  internal  waves  from  the  shipboard  data  collected  on  the 
Endeavor.  The  top  row  shows  cross  sections  at  HH  (left)  and  VV  (right)  polarizations 
while  the  bottom  row  shows  the  corresponding  scatterer  velocities.  Clearly  ao(FiFi)  is 
much  larger  than  ao(VV)  near  the  crests  of  the  internal  waves.  Similarly  scatterer 
velocities  are  much  larger  for  HH  polarization  than  for  W.  Velocities  at  HH  frequently 
exceed  1  m/s,  sometimes  reaching  2  to  3  m/s.  Similarly  high  velocities  can  be  seen  in  the 
VV  image  though  velocities  in  general  are  much  lower  at  this  polarization.  We  interpret 
these  results  as  being  evidence  of  the  involvement  of  breaking  waves  in  microwave 
signatures  of  internal  waves. 
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Figure  12.  Images  of  internal  waves  collected  in  2006  from  the  R/V  Endeavor.  The  top 
row  shows  ao(HH)  (left)  and  ao(VV)  (right);  the  bottom  row  shows  scatterer  velocities  at 
HH  (left)  and  VV  (right).  Normalized  radar  cross  section  scales  are  in  dB;  velocity  scales 
are  in  m/s.  Positive  velocities  are  toward  the  radar.  Dark  blue  in  the  bottom  row 
indicates  that  the  signal-to-noise  ratio  was  too  low  to  determine  a  velocity. 


Unfortunately,  we  were  unable  to  determine  the  phase  of  the  microwave  signatures 
with  respect  to  the  surface  velocities  set  up  by  the  internal  waves  for  these  shipboard 
measurements  because  the  radar  measurements  and  the  ADCP  measurements  did  not 
overlap.  The  radar  operated  early  in  the  experiment  until  a  large  storm  damaged  the 
antennas  and  the  ADCP  did  not  begin  operating  until  after  this  storm. 

One  striking  feature  of  the  images  shown  in  Figure  12  is  that  both  the  cross  sections 
and  velocities  are  highly  variable  within  the  area  of  strong  return.  If  we  average  the  data 
in  these  images  that  are  within  20°  of  the  ship’s  heading,  we  obtain  plots  as  a  function  of 
ground  range  such  as  those  shown  in  Figure  13.  For  the  three  different  times  shown  in 


Figure  13,  the  ship  was  always  heading  into  the  internal  waves  and  the  antennas  were 
always  rotating  eountereloekwise.  Thus  the  VV  antenna  viewed  the  heading  direetion 
about  20  seeonds  later  than  the  HH  antenna  and  the  figure  shows  that  the  internal  wave 
has  advaneed  (moved  left)  during  that  time.  This  motion  ineludes  both  the  IW  speed  and 
the  ship  speed.  All  of  the  panels  of  Figure  13  show  indieations  of  oscillations  of  length 
15  to  30  meters.  For  the  HH  return,  though,  these  oscillations  increase  significantly  on 
the  forward  face  of  the  general  increase  in  cross  section  due  to  the  internal  waves.  The 
internal  wave  increases  the  amplitude  of  these  short  gravity  waves,  leading  to  an 
oscillating  cross  section.  The  differences  between  the  ship’s  heading  and  the  wind 
direction  in  the  three  cases  shown  are  104°,  74°,  and  118°.  Thus,  the  wind  is  almost 
perpendicular  to  the  internal  wave  in  all  cases  so  a  significant  increase  in  short  gravity 
waves  is  expected  in  front  of  the  crest  of  the  internal  wave  as  shown  in  Figure  9d. 

The  increase  in  these  waves  affects  the  HH  cross  section  more  than  the  VV  cross  section. 
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Figure  13.  Average  cross  section  versus  ground  range  for  data  taken  within  20°  of  the 
ship’s  heading.  The  top  row  plots  ao(HH)  as  a  function  of  ground  range  at  various  times 
while  the  bottom  row  shows  similar  plots  of  ao(VV)  at  the  times  that  are  a  few  minutes 
different. 

The  airborne  imagery  also  confirms  the  fact  that  HH  cross  sections  are  more  affected 
by  internal  waves  than  VV.  Figure  14  documents  one  case  of  a  flight  track  across  a  train 
of  internal  waves.  Internal-wave-induced  oscillations  in  the  HH  imagery  are  much  more 
pronounced  than  in  the  VV  imagery  even  though  the  mean  cross  section  level  in  the  VV 
image  exceeds  that  in  the  HH  image.  This  is  also  obvious  in  the  ratio  of  the  standard 
deviation  to  mean  cross  section  also  shown  in  the  figure.  Figure  15  shows  the  same  data 
with  the  mean  cross  section  level  removed  to  emphasize  the  higher  modulation  at  HH 
polarization.  Spectra  of  the  images  in  Figures  15a  and  b  are  shown  in  Figures  15c  and  d. 
We  will  use  below  two  measures  of  both  the  internal  wave  propagation  direction  and  the 
intensity  of  its  surface  signature.  Directions  can  be  obtained  both  from  the  imagery  and 
from  the  spectra.  The  intensity  of  the  surface  signature  is  given  both  by  the  ratio  of 
standard  deviation  to  mean  cross  section  of  Figure  14  and  by  the  rms  level  of  the  image 
spectra  that  are  shown  in  Figures  15c  and  d. 
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Figure  14.  a)  Track  of  the  Skymaster.  b)  Normalized  radar  cross  sections  versus 
incidence  angle  along  the  lines  in  panels  c  and  d.  c)  Image  of  ao(HH).  d)  Image  of 
Oo(VV).  e)  Cuts  through  ao(HH)  image  at  constant  incidence  angles  of  45°,  55°,  and  65°. 
The  numbers  to  the  right  are  mean  cross  section  and  its  standard  deviation  divided  by  the 
mean  at  the  three  incidence  angles,  f)  Same  as  e  but  for  ao(W) . 
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Figure  15.  a)  Zero-mean  image  of  ao(FlH).  b)  Same  as  a  but  for  ao(W).  c)  Spectrum  of 
zero-image  of  ao(HH).  Aircraft  heading  is  shown  above  the  spectrum  and  the  rms  value 
of  the  spectrum,  the  wavelength  of  the  IW,  and  its  direction  of  propagation  in  degT  are 
also  shown,  d)  Same  as  c  but  for  ao(VV).  The  direction  of  the  antenna  look  relative  to 
the  internal  wave  propagation  direction  is  shown  above  the  image. 


On  August  9,  2006,  the  plane  and  the  ship  both  traversed  the  same  internal  wave  train 
many  times.  The  ship  travelled  either  against  the  internal  wave  or  with  it  but  the  plane 
flew  in  many  different  compass  directions.  We  were  therefore  able  to  determine  the 
maximum  NRCS  produced  by  the  internal  wave  looking  both  upwave  and  downwave  at  a 
variety  of  incidence  angles.  From  the  plane  measurements,  we  could  also  determine  the 
dependence  of  the  internal  wave  signature  on  antenna  look  direction.  Figure  16  shows 
the  average  of  the  maximum  NRCS  values  for  both  TUT  and  VV  polarizations  produced 
by  the  internal  waves  when  the  antennas  looked  upwave  (right)  and  downwave  (left). 

Also  shown  in  this  figure  are  the  crosswind  NRCS  values  predicted  by  the  multiscale 
model  (Plant,  2002)  for  a  sea  surface  disturbed  only  by  the  wind.  These  curves  are  the 
same  in  the  left  and  right  columns  since  these  columns  are  the  same  in  the  absence  of 
internal  waves.  We  have  not  shown  multiscale  predictions  for  the  panels  where  the  wind 
speed  is  between  1  and  2  m/s  since  at  these  low  winds,  model  predictions  depend  not 
only  on  the  wind  speed  but  also  on  its  variability.  This  latter  quantity  was  difficult  to 
estimate. 


Figure  16.  Maximum  NRCS  values  of  the  microwave  signature  of  an  internal-wave  train 
as  a  function  of  incidence  angle.  Data  in  the  left  column  were  collected  when  the 


antennas  looked  along  the  internal  wave’s  propagation  direction;  data  in  the  right  column 
were  obtained  with  the  antenna  looking  into  the  internal  wave.  The  lines  are  predictions 
of  the  multiscale  model  of  Plant  (2002)  for  a  sea  surface  disturbed  only  by  wind  averaged 
over  azimuth  angles  between  45°  and  135°.  Black  lines  and  circles  refer  to  W 
polarization;  red  dashed  lines  and  asterisks  refer  to  HH. 

Figure  16  suggests  several  conclusions.  Clearly,  maximum  Oo  values  of  internal  wave 
signatures  are  much  closer  to  multiscale  predictions  at  W  polarization  than  at  HH.  This 
agrees  with  the  less  distinct  internal  wave  signatures  at  VV  than  at  HH  that  we  have 
pointed  out  above.  There  is  some  tendency  for  maximum  ao(VV)  to  trend  closer  to  the 
multiscale  predictions  at  higher  wind  speeds  and  incidence  angles.  This  suggests  that 
internal  wave  signatures  at  VV  polarization  should  be  relatively  difficult  to  detect  under 
these  circumstances.  Maximum  ao(HH)  values  of  internal  wave  signatures  are  much 
higher  than  the  values  predicted  by  the  multiscale  model.  This  model  does  not  include 


wave  breaking  and  we  believe  that  breaking  waves  are  the  most  likely  souree  of  the 
enhanced  backscatter  at  HH  polarization  at  the  maxima  of  the  internal  wave  signatures. 
This  interpretation  is  supported  by  the  observation  that  in  the  shipboard  data,  ao(HH)  is  5 
to  10  dB  higher  than  ao(VV).  We  cannot  rule  out  the  possibility  that  wave  breaking  also 
plays  a  role  in  increasing  ao(VY)  at  signature  maxima  above  its  predicted  value  at  low 
wind  speeds. 

The  flights  of  the  Skymaster  over  all  internal  waves  we  viewed  near  the  shelf  break 
showed  that  there  is  very  little  tendency  for  the  signature  of  the  internal  waves  to  vary 
depending  on  antenna  look  direction  at  moderate  incidence  angles.  Figure  17  documents 
this.  Ratios  of  the  standard  deviation  to  the  mean  cross  section  over  the  internal  wave 
and  the  rms  value  of  the  image  spectrum  are  shown  in  this  figure.  In  all  cases,  no 
dependence  of  the  signature  on  azimuth  angle  can  be  detected.  Note  once  again  that 
internal  wave  signatures  are  more  distinct  in  HH  backscatter  than  in  W. 


Figure  17.  Measurements  of  the  intensity  of  the  internal-wave  microwave  signature  for 
various  antenna  look  directions,  a)  The  ratio  of  the  standard  deviation  to  the  mean  cross 
section  of  the  IW  train  at  an  incidence  angle  of  45”.  b)  Same  as  a  but  for  a  55”  incidence 
angle,  c)  Same  as  a  but  for  a  65”  incidence  angle,  d)  The  rms  value  of  the  image 
spectrum.  This  is  necessarily  calculated  over  the  range  of  incidence  angles  that  formed 
the  image. 

Finally,  in  Figure  18,  we  show  these  same  measures  of  intensity  as  a  function  of  wind 
speed  for  all  the  data  that  we  collected  from  the  plane.  As  is  well  known,  IW  signatures 
become  less  distinct  as  the  wind  speed  increases  and  Figure  18  indicates  this.  We  find 
that  the  signatures  have  nearly  disappeared  at  a  wind  speed  of  9  m/s. 


Figure  18.  Dependence  of  the  intensity  of  internal  wave  signatures  on  wind  speed.  The 
various  panels  show  the  same  dependent  variables  as  Figure  16. 


Conclusions 

We  have  shown  not  only  that  microwave  signatures  of  internal  waves  are  more 
distinct  at  HH  polarization  than  at  VV  but  also  that  maximum  ao(HH)  values  are 
comparable  to  or  even  greater  than  maximum  ao(VV)  values  induced  by  the  internal 
waves.  The  maximum  ao(HH)  values  produced  by  internal  wave  surface  disturbances  are 
much  larger  than  the  multiscale  model  predicts  while  maximum  ao(VV)  values  are 
comparable  to  or  above  the  predicted  values.  Furthermore,  the  intensities  of  the 
signatures  at  both  polarizations  vary  little  with  the  direction  of  antenna  look.  All  of  these 
results  suggest  to  us  that  breaking  waves  play  a  significant  role  in  the  production  of 
microwave  signatures  by  internal  waves.  The  many  ways  that  this  could  happen  were 
outlined  in  Part  1 . 

The  cross  section  maximum  of  the  signature  leads  the  maximum  near-surface  velocity 
of  the  internal  wave  by  larger  distances  the  steeper  the  internal  wave  is.  This  same 
behavior  is  seen  in  simulations  of  the  location  of  the  maximum  modulation  of  short 
gravity  waves  by  internal  waves  using  standard  action  balance  techniques.  This  strongly 
suggests  that  microwave  signatures  of  large  internal  waves  are  produced  by  the  breaking 
of  short  gravity  waves  as  they  are  modulated  by  the  internal  waves.  In  fact,  in  our 
shipboard  measurements  off  the  New  Jersey  coast,  we  observed  cross  section  oscillations 
15  to  30  meters  in  length  that  were  modulated  by  the  internal  waves. 

Most  of  the  microwave  signatures  observed  in  this  low-grazing-angle  study  consisted 
of  increases  of  the  microwave  cross  section  by  the  internal  waves.  A  few  observations  of 
increases  followed  by  decreases  were  also  recorded.  We  did  not  observe  cases  consisting 
only  of  a  decreased  cross  section.  We  note  that  our  simulations  of  short  gravity  wave 
modulation  by  internal  waves  also  showed  cases  of  both  increases  and  increases  followed 
by  decreases  in  the  modulation  intensity.  Da  Silva  et  al.  (1998)  have  reported  seeing 


microwave  surface  signatures  of  internal  waves  from  satellites  at  much  higher  grazing 
angles  that  consisted  of  cross  section  increases,  decreases,  and  paired  increases  and 
decreases.  They  associated  the  first  type  of  signature  with  low-wind  conditions  where  the 
decrease  is  not  visible  due  to  system  noise.  This  was  not  the  case  in  our  low-grazing- 
angle  study.  Our  signal-to-noise  ratio  was  sufficiently  high  that  decreased  cross  sections 
were  not  hidden  by  noise. 

The  movement  of  the  intensity  maximum  down  the  front  face  of  the  internal  soliton 
as  the  wave  steepened  probably  explains  the  results  of  Kropfii  et  al.  (1999)  that  NRCS 
maxima  in  sinusoidal  internal  waves  trains  lie  near  the  minima  of  the  surface  currents 
generated  by  the  IW.  Their  quoted  IW  phase  speeds  are  0.7  to  0.9  m/s  while  the 
maximum  currents  4.4  m  below  the  surface  appear  to  be  less  than  0.7  m/s  and  the  time 
required  for  an  IW  to  pass  a  fixed  point  is  about  0. 1  hours  (full  width  at  half  max).  Thus 
one  thousand  times  the  ratio  of  surface  current  to  width,  a  measure  of  maximum  strain 
rate,  is  about  2.8  1/s  or  a  bit  higher  than  the  highest  found  in  this  study.  Since  wind 
directions  generally  had  a  significant  component  perpendicular  to  the  IW  propagation 
direction  during  the  study,  this  implies  that  their  modulation  pattern  should  lead  the  IW 
current  by  nearly  the  width  of  the  internal  wave  putting  it  nearly  out  of  phase  with  a 
sinusoidal  current.  We  were  not  able  to  verify  this  for  the  sinusoidal  waves  we  observed 
off  the  New  Jersey  coast  because  the  radar  measurements  and  ADCP  measurements  did 
not  overlap. 

Much  additional  information  is  available  in  this  data  set  and  is  being  studied.  In 
particular,  we  are  able  to  obtain  the  speeds  of  the  internal  waves  from  plots  such  as  those 
shown  in  Figures  2  and  3.  We  are  investigating  these  speeds  and  attempting  to  compare 
them  with  predicted  non-linear  wave  speeds.  Finally,  we  have  not  yet  processed  our  data 
to  yield  the  velocity  spread  within  each  range  cell.  This  would  aid  in  the  identification  of 
breaking-wave  events  and  will  be  undertaken. 
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Outline  of  this  paper 


1.  At  low  grazing  angles,  IWs  always  increase  the  cross  section.  (What  about  at  mid 
incidence  angles  from  aircraft?  Can’t  tell.) 

2.  Polarization  ratios  of  IW-disturbed  backscatter  are  often  greater  than  1  for  both 
airborne  and  shipboard  measurements. 

3.  For  shipboard  measurements,  maximum  cross  sections  occur  farther  forward  of  the 
wave  crest  for  steeper  waves,  ie,  farther  forward  on  the  shelf,  near  peak  in  deep  water 

4.  Phase  speeds?  Just  a  little  so  that  location  of  cross  section  maxima  can  be  determined. 
Will  need  IW  widths,  intensities. 

5.  Many  different  kinds  of  surface  signatures  of  IWs  in  the  South  China  Sea.  (Small  ones 
similar  to  those  off  the  NJ  coast?  Next  paper. 

6.  If  possible,  compare  with  theoretical  amplitude  modulations. 

7.  Angular  dependence  of  IW  modulation  from  aircraft  measurements.  (Need  to  get 
ADCP  currents  included) 

Other  Paper: 

1 .  IW  interactions 

2.  IW  phase  speeds  versus  max  velocity 

3.  Compare  with  two-layer  fluid  theory. 

4.  IW  widths,  amplitudes 

Still  Another  Paper  (Like  IGARSSIO  Abstract): 

1 .  Document  environmental  conditions  under  which  breaking  is  important 

2.  Document  radar  conditions  under  which  breaking  is  important 

3.  Include  Doppler  bandwidths  in  the  study. 


